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1. Cutoff effects and renormalization 

Phenomenological results from simulations of lattice QCD to compare with experiments should 
be obtained with all the systematic uncertainties under control. The first requirement is to have an 
efficient algorithm to simulate A'^f = 2 dynamical light quarks with the possibility to include 1(+1) 
heavier quarks. The algorithm should allow, in a reasonable time, to reach small pion masses 
(rrijc < 300 MeV) where a matching with chiral perturbation theory (;^PT) should become pos- 
sible, and to simulate a large enough volume (L > 2 fm). The second requirement is to have a 
lattice action with good scaling and simplified renormalization properties, as close as possible to 
the renormalization of continuum QCD. The topic I want to address in this contribution is if lattice 
twisted mass QCD (tmQCD), combined with a suitable algorithm, is a possible lattice action that 
fulfills these requirements. 

1.1 Lattice QCD action 

Despite the only rather recent interest, the tmQCD fermionic lattice action has a long history. 
It was introduced in [[I]] as a tool to study spontaneous parity and flavour symmetry breaking. In 
it was proved that lattice tmQCD is an alternative discretization of lattice QCD. The lattice QCD 
action 

S = Sf,[U]+SF[U,\lf,W] (1-1) 
has a fermionic part given by tmQCD 

=«^L{¥W[£»[t/] +'«o + 'M75T^]'K(x)} (1.2) 



and for the moment we leave unspecified the gauge part Sg. In eq. ( [L^ ) D[U] is the massless 
Wilson-Dirac operator 

D[U] = l[YA^n + y;)-aV;V^,] (1.3) 

/Mo is the untwisted bare quark mass parameter, /i is the bare twisted quark mass, and is the 
third Pauli matrix acting in flavour space. In it was shown that the standard framework of 
the Symanzik improvement program, works in the similar way as for usual Wilson fermions. In 
particular for spectral quantities no further improvement coefficients are needed. A set of scaling 
tests have been performed, using the non-perturbatively improved clover action with twisted mass, 
in small [Q] and large [Q] volume, confirming that the usual Symanzik improvement program can 
be applied also for tmQCD. 

In a remarkable paper of Frezzotti and Rossi a step forward was made. It was proved that 
parity even correlators of multiplicatively renormalizable fields, are free from O(fl') effects, and so 
no improvement coefficients are needed (automatic 0(a) improvement), if the target continuum 
theory is fully twisted ^ . The proof in is based on a set of spurionic symmetries of the lattice 
action. Here we give a simpler proof based on the symmetries of the continuum QCD action (see 
appendix A of [Q] for an analogous proof). The Symanzik ^ ^ [T^, |ll|, 12] effective action reads 



>eff ■ 



■.So + aSi + ... (1.4) 



' To obtain automatic 0(a) improvement in ||^ also other possibilities were exploited. Here we will concentrate on 
the automatic 0(a) improvement that is used in numerical simulations. 
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and we are interested in a continuum target theory where the physical quark mass is fully given by 
the renormalized twisted mass /Ir (fully twisted theory) 

So = J d'^xxj/ix) [y^D^ + iAiR/s^^] vi^) (1.5) 

The correction terms in the effective action are given by 

S, = / d^y^, (y) ^1 (y) = Y^c^y) (1.6) 

•' i 

where the dimension five operators classified on the basis of the symmetries of the lattice action 
are given by 

ffi = \j7a^yF^yY ^2=At^W ^3=A^W (1-7) 

where A is an energy scale of the order of the QCD scale Aqcd- The operator ffi is the usual clover 
term. The operators ^2 and 1^3 are related to the renormalization of the untwisted quark mass. 
Since we are interested in a continuum target theory where the untwisted quark mass vanishes, 
the operator 1^3 parameterizes the mass independent 0(a) uncertainties in the critical mass. We 
consider now a general multiplicatively renormalizable multilocal field that in the effective theory 
is represented by the effective field 

^efi = ^o + a^i + ■■ (1-8) 
A lattice correlation function of the field to order a is given by 

(<!>) = {^o)o + a I J^(<Doifi(j))o + a(<J>i)o + ... (1-9) 

where the expectation values on the r.h.s are to be taken in the continuum theory with action Sq. The 
key point is that the continuum action (jTsb is symmetric under the following parity transformation 



Y{x) — > ro{ir5T:^)¥{xo, -x) (1.10) 

\j7{x) — >\j7{xo,-x){iY5T^)Yo (1.11) 

and that all the operators in eq. ([TtI), of the Symanzik expansion of the lattice action, are odd under 
the parity symmetry of the continuum action. If the operator <t>o is parity even, the second term 



in the r.h.s. of eq. (1.9) vanishes, and <I>i, being of one dimension higher, is parity odd: for the 
same reason the third term in the r.h.s of eq. (|L^) vanishes. Possible contact terms coming from 
the second term amount to a redefinition of and so do not harm the proof. 

It is then also clear that in order to achieve automatic 0(fl;) improvement, the continuum target 
theory must have a vanishing untwisted quark mass hzr, otherwise the standard mass term 
will break the parity symmetry of the continuum action defined before. The most natural way to 
achieve this on the lattice is by setting the untwisted bare quark mass to its critical value niQ = mc. 
The proof also shows that a possible uncertainty of 0(a) in the critical mass does not wash out 
automatic 0(a) improvement since these uncertainties, are odd under parity. A remark is in order 
now. We take the polar mass defined in ^ 



M = ^1 ^'^ + m\ = \l jX^ + {r]iaK^Y; m^ = mQ-mc (1.12) 
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where the 77 1 term parameterizes the mass independent 0(a) uncertainties in the value of the un- 
twisted quark mass m^. Expanding in powers of a we have 



(1.13) 



We observe immediately that as soon n < aA^, even if parametrically 0(a) terms are absent in 
( 1.13 ), there is a term of O(a^) with a coefficient that tends to diverge as soon /i is made smaller and 



smaller. From this example we can conclude that to have an effective automatic 0(a) improvement, 
without big O(a^) effects, with a generic choice of the critical mass, such that the uncertainties 
in the untwisted quark mass are of order aA^, we need to have the constraint n > aA^. It has 
been shown in [^] that these cutoff effects that diverges at small quark masses, so called infrared 
divergent (IR) cutoff effects, are a general property of tmQCD. These dangerous cutoff effects are 
removed by an appropriate choice of the critical mass. 

1.2 0(a) improvement and small pion masses 

The 0(a) uncertainties of the untwisted quark mass depend on how the critical line is fixed, 
hence the choice of the critical mass has to be discussed with care. The issue was raised by the 
work of Aoki and Bar [|lj] and by the numerical results obtained in [14]. This problem has been 
further analyzed in several aspects [jTs], 0, |l6||. In [13, 15, 16] the theoretical framework is twisted 
mass chiral perturbation theory (tm;^PT) [jlTI] where the cutoff effects are included in the chiral 



lagrangian along the lines of [[18|, [19|]. In this framework a power counting scheme that includes 
quark mass and lattice spacing has to be specified. In particular in [ |T3| ] the power counting was 
/i ~ a^A^ while in [15] it was jj. ~ aA^. We stress here that this approach for the description of 
lattice data, does not require a continuum extrapolation, hence the power counting scheme does not 
mean that /i goes to zero in the continuum limit but represents only an order of magnitude equality. 



Both these works [13, 15] agree on the fact that choosing the critical mass imposing a vanishing 
PCAC quark mass 

I,(5oAg(x)P«(0)) 



"JPCAC 



2Ix(^''W^"(0)) 



1,2 



where 



(1.14) 



(1.15) 



(1.16) 



allows to have automatic 0(a) improvement, and in particular down to quark masses that fulfill 
fX ~ a^A^ for [ [l3| ] and a^A^ < jJ. < aA^ for [15]^. In [Q] a Symanzik expansion along the lines of 
[12] was performed confirming the results of [13, [I^ ]. 

A possible practical procedure is then to compute for a fixed value of jj. the critical mass 
rric from the vanishing PCAC mass, and then to extrapolate the set of critical masses obtained 
for different values of /i to /i = (method A). This procedure has been used in [|2^, 21]. In 
fig. |l| a typical extrapolation of the critical mass to = is shown. With this procedure the 0(a) 

^We will see in section Isii] that the phase structure of N{ = 2 dynamical Wilson fermions does not allow anyway the 
twisted mass to be smaller then fj^ ~ cP'A^. 



014/4 



Twisted mass lattice QCD 



Andrea Shindler 



uncertainties of the critical mass are fixed in such a way that, for a generic value of jj., the dangerous 
aA^ cutoff effects in the untwisted quark mass are absent. The slope of the curve is proportional, as 




Figure 1: Determination of the critical mass nic (fc^r' = lanic + 8) for given values of at j3 = 6.0, and 
extrapolation to /i = 0. The red point is the critical mass determined with method C (see text). The difference 
between the two determination of the critical mass should be an 0{a). 



it has been recently discussed in [23, 16], to 0(a) cutoff effects related to the discretization errors 



of the PCAC mass. We remind that it is not surprising that the PCAC mass is not automatically 
0(fl;) improved since it is an odd quantity under the parity transformation of eqs. ( |1.1(]| , |1.11[ ). 

In [ 23 1 the extrapolation to = is not performed and each value of the critical mass has been 
used for the corresponding value of /i used in the simulations (method B). With this method the 
0(a) cutoff effects of the critical mass are obviously fixed in such a way that the untwisted quark 
mass is always vanishing for all the simulation points. 

These two methods, even if they give different cutoff effects to the critical mass, are perfectly 
good in order to achieve automatic 0(a) improvement. Another possible way to fix the critical 
mass, expecially practical for expensive dynamical simulations, is to compute the critical mass, 
using the PCAC relation at the smallest value of /i, and then use this critical mass for all the 
simulation points at heavier masses. 

Using methods A and B a set of quenched studies [24, 14, 25, 21, ^ have been 

performed to check the result of ^ and to gain experience with this formulation of lattice QCD. 
An interesting quantity to compute with tmQCD is the pseudoscalar decay constant /ps. As it 
was noted in ||5|, 24], the computation of /ps does not require any renormalization constant, in 
contrast of ordinary Wilson fermions, and moreover given automatic 0(a) improvement, does not 
need the computation of any improvement coefficient. Thus the situation for this quantity is like 
with overlap fermions. In fig. |^ (left panel) the continuum limit of ro/ps ^, the critical mass being 
computed with method A, is shown as a function of (a/ro)^. The scaling is consistent with being of 
O(a^), and moreover the O(a^) effects are rather small for all the pseudoscalar masses investigated 



'The values of ro/a, f-Q = 0.5 fm being the Sommer scale [E9||, are taken from |pO|| 



014/5 



Twisted mass lattice QCD 



Andrea Shindler 



■ 


K "^""^ 




=272 MeV 


• 


K '^"'^ 




=514 Me V 


A 


^ PCAL 




=900 MeV 


1,1, 




0.05 




0.06 



(a/r„) 



0.18 
0.16 
0.14 



T — ' — I — ' — r 

: fps[GeV] 



-| — ' — I — ' — I — ' — I — ' — r 




• Wilson imQCD at Ji/2 
□ NP 0(a) improved Wilson 



nip/ [GeV 



0.2 0.4 0.6 0.8 



Figure 2: Left panel: scaling behaviour of ro/ps for 3 fixed values of romps. Right panel: Continuum limit 
values for /ps as a function of nip^ in physical units. The empty squares are taken from [^. 



down to mps = 272 MeV. The right panel of fig. ^ shows the chiral behaviour of the continuum 
pseudoscalar decay constant, compared with the non-perturbatively 0(a) improved data of [p7|]. 
We remark that this comparison is purely illustrative since it is in the quenched approximation, and 
the simulations with clover fermions had to stop around mps ~ 500 MeV due to the appearance of 
exceptional configuration. To see the potential of tmQCD another interesting phenomenological 
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Figure 3: {x)^^{lJ. — 2 GeV) extrapolated to the continuum as a function of the pion mass. Open squares 
represent results that are obtained from a combined continuum extrapolation of earlier Wilson and clover- 
Wilson simulations [|T|. The filled circles represent results using Wilson twisted mass fermions [p2|]. The 
open circle denotes a result which is not corrected for finite size effects and the diamond corresponds to the 
experimental point. 



quantity is the average momentum carried by valence quarks in a pion ((x)). In [32] results using 
tmQCD were presented. Here we concentrate on the chiral behaviour in the continuum, having in 
mind that the renormalization has been performed already in a non-perturbative way [51, 31]. Fig. 



014/6 



Twisted mass lattice QCD 



Andrea Shindler 



|3| shows that in principle also for this quantity small pseudoscalar masses mps < 300 MeV can be 
reached, opening the possibility of a safe chiral extrapolation. 



P=6.0 
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0.05 0. 1 



0.2 



0.25 



0.3 



(ma) 




Figure 4: Left panel: comparison of the chiral behaviour at fixed lattice spacing (j3 = 6.0) of the pseu- 
doscalar decay constant computed using method A, B, C and with results obtained with overlap fermions. 
Right panel: unconstrained continuum limit, for several values of fixed charge pion masses, of ro/ps per- 
formed using method A and C to determine the critical mass. 

In [U to obtain automatic 0(a) improvement the critical mass nic was computed extrapolating 
the squared pseudoscalar mass to the chiral limit using data from the pure Wilson theory (method 
C). Using this determination of nic several quantities were computed. In particular in the left 
panel of fig. ^ there is a comparison of the chiral behaviour at fixed lattice spacing (j3 = 6.0) of 
the pseudoscalar decay constant computed using method A, B, C and with results obtained with 



overlap fermions [|34|]. While methods A, B and the overlap data are all consistent within the 
statistical errors,"^ the data obtained using method C to fix the critical mass, show a "bending" 



towards the chiral limit. The same phenomenon was observed also for the vector mass []14|]. The 
"bending" phenomenon appeared exactly when /i ~ aA^. Having in mind the caveat observed 
before in the proof of automatic 0(a) improvement, this indicates that the extraction of the critical 
mass with method C leaves the dangerous aA^ in the untwisted quark mass uncanceled. This is 



numerically confirmed by the results of [21 1, showed in the right panel of fig. ^, since using method 
A and C to determine the critical mass, a consistent continuum limit is obtained, showing also that 
method C induces big O(a^) effects and a reduced scaling window. 



A description of the "bending" phenomenon at fixed lattice spacing, has been obtained in [22] 
using ;tPT, as it is shown in the left panel of fig. B|, where a fit to available quenched data is per- 

2 2 

formed on the ratio R = — This analysis shows also that ^PT theory is able to describe the 



lattice data up to ~ 80 MeV. It is reassuring that using method A to determine the critical mass 
and restricting the data to the region were ;^PT is applicable the ratio R is flat (right panel of fig. 
^ consistently with continuum ;^PT (up to chiral logs). In [^], based on the observation that the 
big O(a^) effects come from uncanceled 0(a) of the PCAC mass, to eliminate the "bending" phe- 
nomenon has been proposed to use a non-perturbatively improved tmQCD action. This approach 



We recall here that since the comparison is made at fixed lattice spacing the data in principle could disagree due to 
different cutoff effects. 
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Figure 5: Left panel: bending phenomenon at j3=6.0 on the ratio R 
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and its description with ;ifPT. 



Right panel: comparison of the ration R using method A and C to determine the critical mass. 



has been numerically tested in [35], and as it can be seen in fig. |6| indeed it confirms that the 
"bending" phenomenon also in this case it is not present. 
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Figure 6: Comparison of the chiral behaviour at fixed lattice spacing (j3 = 6.0) of the pseudoscalar decay 
constant computed using method A and C for both tmQCD and non-perturbatively improved tmQCD. 



1.3 Renormalization 



In [36] it has been given a construction of a Schrodinger functional (SF) with twisted bound- 



ary conditions that preserves the nice properties of 0(a) improvement without bulk improvement 



coefficients; (see [37] for a possible alternative to this construction). The construction is based 



on the consideration that in a finite volume with suitable boundary condition the Wilson theory 



in the chiral limit is 0(a) improved, and it makes use of orbifolding techniques (see [ ]38[ ] for an 
application of orbifolding techniques to Ginsparg-Wilson fermions). 

A simple way to visualize the construction is to repeat the proof of automatic 0(a) improve- 
ment given in section |1.1| , where now since we are in a finite volume with suitable boundary con- 



ditions, the twisted mass could be safely sent to zero. Then the new boundary projectors [36] 
Q± = ^(1 +170 75 T^) commute with the previous parity transformation (as the twisted mass term 
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in infinite volume). It is very important to note that the new boundary projectors can be obtained 
performing a chiral rotation of the original projectors in the standard SF framework [|^]. An im- 
portant consequence is that the running of the coupling constant, should be identical to the running 



computed with the "old" SF [|40|]. The 0(a) uncertainties in the critical mass do not harm the 0(a) 
improvement. 

2. Flavour symmetry 

When tmQCD is used to define the standard QCD correlation functions some of the physical 
symmetries are restored only in the continuum limit. In particular flavour and parity symmetries. 
The explicit breaking of flavour symmetry generates for example splitting between charged and 
neutral pions, while the absence of parity symmetry, gives as a consequence the appearance of 
states of opposite parity in the spectral decomposition of usual correlators. Both these phenomena 
are expected to vanish, at maximal twist, with a rate of O(a^) pTl]. Here we concentrate on the 
flavour symmetry breaking. 

To fix the notation we recall some basic definitions. The charged pseudoscalar currents are 
given by 

P^{x)=W{x)y, — ^if{x) T± = ^— (2.1) 
and a possible interpolating field for the neutral pion is the scalar current 

^{x) = ^^{x)\if{x). (2.2) 

The charged and neutral pseudoscalar masses can be extracted by the following correlators 

Q+(xo)=a^^(/'+(x)p-(0)) Qo(xo) = a^£(S«(x)50(0)) (2.3) 

X X 

Qo(xo) =a^£{(-tr[G(0,x)G(x,0)] +tr[G(x,x)]tr[G(0,0)])} (2.4) 



where G(x,j) is the fermionic propagator. In [ ]42[ ] a pilot quenched study has been preformed to 



study flavour breaking effects with tmQCD. For the neutral pseudoscalar correlator in eq. (2.4) a 
first possibility is to study only the connected part. In the quenched approximation it is still possible 
to interprete the connected part in terms of local operators. The reason is that one could think the 



connected part as coming from Wick contractions obtained using the Osterwalder-Seiler (OS) [43] 
action 

Sos =a'Y. {W{x) [D[U] + mo + ipLJiWi^) } • (2.5) 

X 

This action has a trivial flavour structure and so does not present any flavour breaking, and in partic- 



ular the disconnected part of eq. ( [2.4] ) vanishes. We remark that this is not the neutral pseudoscalar 
meson of tmQCD, but it is an interesting quantity to study with precise data on its own, in view 
of a possible use of mixed actions (the OS action for the valence quarks and tmQCD for the sea 
quarks). 

In fig. ^ the scaling behaviour of the connected correlator (OS pseudoscalar), is compared 
with the pseudoscalar meson in tmQCD where also the disconnected part is included (in both 
computations method A is used for the determination of the critical mass). For all the technical 
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0,08 



Figure 7: Scaling behaviour of the mass splittings between the neutral and the charged pseudoscalar masses 
for 2 values of romps- The open squared are the data for the neutral pseudoscalar meson with tmQCD, and 
the stars only the connected contribution (pseudoscalar meson with OS action). The full and dotted lines are 
an estimate of the a^ dependence for the two pion splittings, making the hypothesis that O(a^) effects are 
mass independent. 



details of the computation of the disconnected part I refer to 44]. The results show an 0{a^) 
scaling for both the pseudoscalar masses, even if there are indications that the neutral pseudoscalar 
meson for tmQCD (with the inclusion of the disconnected correlator) has reduced cutoff effects, 
within the rather large statistical errors. It is possible to give a very rough estimate of the pion 
splitting tf^{m^o ~^\±) — c{a/rQY with c ~ 10 (with large errors). Comparing to a quenched 



simulation for naive staggered fermions with Wilson gauge action [45], one finds a similar size of 



the flavour splitting encountered for the pion mass at a similar lattice spacing with a value c ~ 40. 



For dynamical improved staggered fermions a value of c ~ 10 has been found [ ]46| ] 



An interesting study of the flavour breaking effects was presented at this conference in []260. To 



avoid the computation of disconnected diagrams in the quenched approximation a second doublet 



for strange and charm quarks is introduced following the strategy of ||47|]. Then the splitting on the 
kaon system is studied. In this study method B has been used for the determination of the critical 
mass. The results shown in fig. indicates that, as expected, the flavour breaking effects vanish 
linearly with a rate of a^, but that indeed they could be significant at a lattice spacing a > O.lfm. 
Recent results with Nf = 2 dynamical tmQCD fermions and DBW2 gauge action (see next 
sections for details on the simulation parameters), indicate that at a lattice spacing of a ~ 0.12 fm 
and a mass /i ~ 12 MeV, the pion mass splitting even if with large errors, is consistent with zero. 

3. A^f = 2 

3.1 Algorithmic improvements 

In the Lattice 2001 conference A. Ukawa presented [SI] a rather impressing analysis on the 
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200 




Figure 8: Scaling behaviour of the mass splitting between neutral and charged kaons. 



possibility of simulating light quark masses with Wilson fermions. This was summarized with 



the now well known Berlin wall figure (see [ |49| ] for a recent update). Recently new algorithms 



[ pOj plU have been proposed that have finally moved the wall to rather small quark masses. Both 
the algorithms are based on the standard HMC but have used new preconditioned In [|^] it was 
shown that with a domain decomposition (DD) preconditioning combined with a multiple time (mt) 
scale integrator [52], light quark masses (rriji = 294 MeV) are reachable with Wilson fermions with 



remarkable performances. In [51] another very efficient preconditioner for the HMC algorithm 



has been introduced and tested, based on a mass preconditioner [ ]53[ ] (also known as Hasenbusch 
(H) acceleration) with again a multiple time scale integrator. In table [I] is summarized a rough 
comparison between the two algorithms, using different lattice actions, based on the so called 
cost figure v = 10^^ {2N + 3)rint{P) introduced in [pO[l. The conclusion is that the algorithms 
have comparable performance down to pion masses of the order of m^i ~ 300 MeV. In fig. ^ is 
plotted the update of the Berlin wall figure. On the left panel there is a comparison between the 
results of [^, 54] (squares and diamond) and the results of [^] (circles). The lines are functions 
proportional to (mps/my)"^ (dashed) and (mps/my)^ (solid). On the right panel it is shown a 



comparison between the Ukawa's formula in [ |48| ] (solid line) and the extrapolation of the results 
in [51] using a (mps/my)^ (dashed) and a (mps/my)^ (dotted) dependence for the data. The arrow 
indicates the physical pion to rho meson mass ratio. In addition there are also data points from 
staggered simulations (see [ p9| ] and references therein). In particular this plot indicates that running 
for one year a 1 Tflop sustained performance machine allows to generate at the physical point with 
a ~ 0.08 fm and a lattice of 24^ x 40, 1000 independent trajectories. 



Action 


Algorithms 


ro/a 


niji [MeV] 


V 


■^int 


W+W 
tlSym+Wtm 


(mt)(DD)HMC [50] 
(mt)(H)HMC [^] 


6.40(15) 
5.20(25) 


294 
280 


0.74(18) 
0.49(34) 


21(5) 
21(14) 



Table 1: Comparison of the 2 algorithms discussed in the text for a similar physical situation. 
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Figure 9: Computer resources needed to generate 1000 independent configurations of size 24^ x 40 at a 
lattice spacing of about 0.08 fm in units of Tflops • years as a function of mps/m\/. See text for a detailed 
description of the plots. 



3.2 Phase diagram of Wilson fermions 



o 
o 

( n 



In [56] the first study of tmQCD with Nf = 2 dynamical fermions was performed. Starting 
the exploration of a completely new territory, it is always good to remember a sentence of G. 
Parisi [^] "Let me describe a typical computer simulation: the first thing to do is to look for 
phase transitions". It is important to have then the correct understanding of the phase diagram with 
Wilson fermions in the 3 parameters space (j3 = 6/go)"^o,M)- To check that the results are not 
induced by the algorithm used it is always good to have at least 2 algorithms that reproduce the 
same results. 

Indeed in [56] using the so called TSMB [58] and GHMC [|^, algorithms rather surprising 
results were found. The action used was Wilson gauge action combined with Wilson fermions with 
and without twisted mass. In particular at a lattice spacing of a ^ 0.16 fm, strong evidence of a 
first order phase transition was found for a rather large range of values of twisted masses going 
from zero twisted mass to /i ~ 100 MeV. This study reveals also that the phase transition tends 
to disappear increasing the value of jj., it persists for = and it is volume independent. A 
typical example of a MC history for the plaquette expectation value can be seen in fig. [l^, where 
a cold and a hot start was performed. These results can help to see from a different point of 



view old numerical and theoretical works. In [|60|, |61j] from a finite temperature study there was 
an indication of difficulties in observing a phase with spontaneous breaking of flavour and parity 
symmetry (Aoki phase) at /3 > 4.8. In [ ]62| ] the MILC collaboration found a surprising bulk first 
order phase transition for Wilson fermions at j8 ~ 4.8. In [ ]63| ] an analysis using the linear sigma- 
model is performed, finding an indication of two possible patterns of symmetry breaking at finite 



lattice spacing. This observation was put on firmer theoretical basis in [18]. In this very important 
paper several interesting results and consideration were done, that, seen now from a different point 
a view, can help to understand the rather surprising numerical results obtained in ]j5^]. In [|T^] for 
the first time the concept of chiral lagrangian at finite lattice spacing is given. The key point of the 
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50 100 150 200 250 300 350 400 450 500 
*HMC 

Figure 10: Metastable states at j3 =5.2. Left panel: MC history of the average plaquette value with a twisted 
quark mass ~ 10 MeV and a lattice size 16^ x 32. Right panel: MC history of the average plaquette for 
pure Wilson fermions (ji = 0) and a lattice size 12^ X 24. 

O 

construction is the observation that the Pauli term in the effective Symanzik lagrangian transforms 
under chiral rotation exactly as does the mass term. I would like to add that this is also the key 
point for the automatic 0(a) improvement for tmQCD at maximal twist. Neglecting the derivative 
interaction, being interested in the vacuum state, the potential of the effective chiral lagrangian 
reads 



4 ^ ' ' 



C2 ~ m'^K^ + m'aK^ + a^K^ 



m 



m — ah 



(3.1) 
(3.2) 



where £ is the matrix that collects the Goldstone boson fields of the theory. We remark here that 
m' is a redefinition of the untwisted quark mass that includes the 0(a) coming from the clover 
term. Up to O(a^) m' is proportional to the PCAC quark mass. The two terms in the potential 
become comparable when m' ~ a^A^. In this region of quark masses the competition of these two 



terms causes a non-trivial vacuum structure that gives the following 2 scenarios [ ]18[ ] for the phase 
diagram of Wilson fermions: 1) The Aoki phase [[l]]; 2) The existence of a 1**' order phase transition 



[18]. The extension to tmQCD of these results is done in [64, 65, 66]. The result is summarized 



in fig. ^ where the x-axis is m'/a^ and the y-axis is /i/a^. A non-zero value of the twisted mass 
washes out the Aoki phase introducing an explicit breaking of flavour and parity symmetry. As can 
be seen from fig. |n| (left panel) the Aoki phase lies on the untwisted axis. In the second scenario 



in fig. [11| (right panel) the first order phase transition line extends into the twisted direction to a 
distance of /Xc ~ a^ A^ . The transition ends with a second order phase transition point, where the 
neutral pion mass vanishes. Several comments are in order now. The occurrence of one of the 
two scenarios depends on the sign of the coefficient C2 proportional to the O(a^) term in the chiral 
lagrangian. This coefficient C2 depends on the choice of the gauge action, on the presence in the 
lattice action of the clover term and on the bare gauge coupling. An analysis with Wilson fermions 
of the two dimensional Gross-Neveau model [|^] indicates that indeed both the scenarios describe 
the phase structure of Wilson fermions depending on the value of the couplings of the model. The 
analysis shows that at strong coupling there is an Aoki phase while at weak coupling the first order 
phase transition line sets in. This analysis has been recently extended for the twisted mass case 
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Figure 11: Left panel: the phase diagram of Wilson diagram according to X PT for C2 > 0. Right panel: as 
the left panel but for C2 < 0. The x-axis is m' / a^ and the y-axis is jx /a^. 



58p, indicating even more complicated structures, like a coexistence of tlie two scenarios at the 
same value of the coupling. 

Our present understanding of the lattice QCD phase diagram can be summarized as follow- 
ing. For values of the lattice spacing much coarser than a = 0.15 fm there is a second order phase 



transition from the standard lattice QCD phase to the Aoki phase |[]J, |69|,|70|]. For smaller values of 
the lattice spacing a first order phase transition appears [|^, ^ ^] that separates the positive 
quark mass from the negative quark mass phase. This first order phase transition is reminiscent 
of the continuum phase transition when the quark mass is changed from positive to negative val- 
ues with the corresponding jump of the scalar condensate as the order parameter of spontaneous 



chiral symmetry breaking. The generic phase structure of lattice QCD is illustrated in fig. 12 and 



discussed inrefs. [pq, 71, 72] 




I Aoki phase 

I 1^ order phase transition plane 



P 



Figure 12: Current knowledge of the Wilson lattice QCD phase diagram as function of the inverse gauge 
coupling j3 o= 1/^^, the hopping parameter if and the twisted mass parameter jX. 



3.3 Minimal pion mass 

In the scenario with a first order phase transition the pseudoscalar mass mps cannot be made 
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j8 =3.65 


p = 3.75 


j8 = 3.90 


a^t =0.01 
a ^0.13 fm 
Lr^ 1.56fm 
K)min~450 MeV 


ajj. = 0.0094 - 0.005 

aR^0.12fm 

L«2fm 

K)mm~400MeV 


ajJ. = 0.0075 - 0.004 
a^iOA fm 
L?a 1.6fm 

{m^)min-2m MeV 



Table 2: Summary of the simulation parameters for dynamical runs of tmQCD with tlSym gauge action. 
The last line is an estimate of the minimal pion mass reachable without encountering metastabilities. 



arbitrarily small both if the chiral point is reached from the untwisted or the twisted direction. 
Lowering the quark mass from the untwisted direction the algorithm will start to sample also in the 
region with negative masses. The minimal pion mass reachable will then depend on the algorithm 
used and on the strength of the phase transition. Lowering the quark mass from the twisted direction 
there is a minimal pion mass given directly by the extension of the first order phase transition line, 
even if the twisted mass gives a sharp infrared cutoff in the sampling performed by the algorithm. 

It therefore becomes important to understand the phase structure of lattice QCD as a pre- 
requisite before starting large scale simulations. As we have seen the extension of the first order 
phase transition line in the twisted direction is proportional to the coefficient [c2|. This coefficient 
depends both on the gauge action used and on the presence of the clover term in the lattice action. 
In [ |73| ] has been studied the lattice spacing dependence of of the first order phase transition with 
Wilson gauge action, taking as a measure of its strength, the gap between the two phases in the 
plaquette expectation value and in the PCAC quark mass. The qualitative estimate for the lattice 
spacing, where a minimal pion mass nij^ ^ 300 MeV could be reached, without being affected by 
the first order first transition is 0.07-0. 1 fm. 

It is suggestive that at the microscopic level the occurrence of this first order phase transition is 
accompanied by a massive rearrangement of the small eigenvalues of the Wilson-Dirac. This rear- 
rangement could be suppressed by the use of a renormalization group improved or O(a^) improved 



gauge actions, and indeed results from [ ]74| ] indicate that metastabilities in the average plaquette 
observed for Nf = 3 dynamical Wilson fermions with a clover term (there is also an indication that 
the same metastabilities survive without a clover term for A'f = 3), can be suppressed replacing the 



Wilson gauge action with the Iwasaki action [75]. 



3.4 Tree-level Symanzik improved gauge action 

The dependence of the phase diagram on the gauge action used and on the lattice spacing 



has been studied in a set of papers [pq, 71, 72, [73n (see also nJa] for a detailed summary of these 



results). The gauge actions so far studied can be parameterized by 

SG = P[bo I (l-V><i(x;At,v))+Z,i(l-Vx2(^.^^^)^] ^33) 

with the normalization condition bo = I — Sbi. The parameters of the tree-level Symanzik action 
[ ]77[ ] (bi = —J2) simulations are summarized in tab. |2[ The last line indicates an estimate of the 
minimal pion mass reachable at the corresponding lattice spacings. In order to check for a possible 
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phase transition and corresponding metastabilities a measure of the average plaquette value as a 
function of the hopping parameter K on runs that start from both a hot and a cold configuration 
has to be done. Since the metastability, if any, will show up around jq. (determined monitoring the 
PC AC mass mpcAC at the corresponding fixed value of /x) attention should be given to the hot and 
cold runs on Kr-values closest to Kc only. 

At j3 = 3.65, a ~ 0.13/m, 12^ x 24 and /X ~ 15 MeV there are signs of a very nearby phase 
transition, as can be deduced from the steep rise in K of the plaquette expectation value (left panel 



in fig. 13 ), from a very slow thermalization and large fluctuations of the plaquette MC history value 



over several hundreds of trajectories (right panel in fig. 13). An estimate of the pseudoscalar mass, 
close to Kc,is mps — 450MeV. 




1000 1500 
MC time 



Figure 13: Left panel: Average plaquette value vs. K" at j8 = 3.65, ro/X = 0.038 on a 12^ x 24 lattice from hot 
(red symbols) and cold starts (blue symbols). Right panel: Average plaquette MC time history for two runs at 
j3 = 3.65, ro^U = 0.038, k = 0.17024 on a 12^ x 24 lattice starting from hot (red line) and cold configuration 
(blue line). 



At j3 = 3.75, a ~ 0.12 fm, 12^* x 24 and ~ 8 MeV there is a similar situation we have 
observed before at j8 = 3.65 and ~ 15 MeV. This is described by fig. (left panel), where it is 
plotted the k dependence of the PCAC mass. This dependence is very useful to monitor a possible 
metastable critical point, since this shows up in a different extrapolated Kc, when the extrapolation 
is performed from positive or negative quark masses. A second twisted mass ~ 15 MeV has 
been simulated in a lattice 16^ x 32 around the critical point for this lattice spacing. Even if a 
strict check done with a hot and a cold start is not available at the moment the K dependence of the 
PCAC mass for this second value of n suggests that the critical point is free from metastabilities. 
The pseudoscalar mass measured for the heaviest twisted mass is around mps ~ 400 MeV. At 
j3 = 3.9, a ~ 0.1 fm, 16^ x 32 and /i ~ 8 and 15 MeV, there are no signs of metastabilities at the 
two corresponding critical points. In fig. |l^ (right panel) is plotted the K dependence of the PCAC 
quark mass. The pseudoscalar masses obtained for the two values of n are respectively mps ~ 280 
and 450 MeV. We remark also that the physical volume at this j3 value is rather small L ~ 1.6 fm, 
and results obtained for pure Wilson fermions [50, 55], indicate that for these quark masses and 
these volumes the finite size effects could be substantial. The estimate of the minimal pseudoscalar 
mass for this lattice spacing is then clearly only an upper bound. 
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Figure 14: PCAC quark mass mpcAC vs. )f on a 16^ x 32 lattice. Left panel: a ~ 0. 12 fm, ~ 8 MeV. Right 
panel: a ~ 0.1 fm, /i ~ 15 MeV. 



15 = 0.67 


j8 = 0.74 


a^u =0.01 
a?s0.19fm 
L 2.3 fm 
Mmin ~ 360 MeV 


an = 0.0075 
a ^0.12 fm 

{m^Uin ~ 320 MeV 



Table 3: Summary of the simulation parameters for dynamical runs of tmQCD with DBW2 gauge action. 
The last line is an estimate of the minimal pion mass reachable without encountering metastabilities. 



3.5 DBW2 gauge action 

In this section I summarize the results obtained using the so called DBW2 gauge action 
[78] {bi = —1.4088). The parameters used in the simulations are summarized in tab. ^ The 



twisted mass for the two lattice spacing is kept roughly fixed to jLt « 12 MeV. The last line indicates 
an estimate of the minimal pion mass reachable at the corresponding lattice spacings. Also for this 
gauge action several quantities have been computed. Here we concentrate as before on the PCAC 
mass and on the minimal pion mass. In contrast with the tlSym results here simulations at full twist 
were never performed, so the evidence for a metastability region can be deduced only indirectly 



from the dependence of the PCAC mass on the untwisted quark mass as discussed before. In fig. |15 
is shown the 1 / {2k) dependence of the PCAC mass for the two lattice spacing used. At a « 0. 19 
fm there is an indirect evidence of a small metastability at full twist, that seems to disappear at 
a Rr: 0.12 fm. To summarize, there are first indications, that in order to reach pion masses of the 
order of mps ~ 300 MeV with tmQCD a gauge action like tlSym or DBW2 is appropriate. There 
are also first indications that this pion mass can be reached with DBW2 at slightly coarser lattices. 



In order to avoid possible large cutoff effects with DBW2 (see for example fig. 3 in |7S]), or 
big coefficients in perturbative expansions, with the present data, the tlSym gauge action may be 
considered a better choice. 

3.6 A^f = 2 + 1 + 1 

The fact that tmQCD can be formulated only for an even number of flavours is not a limitation. 
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Figure 15: Determination of the critical hopping parameter Kq by extrapolating to zero the untwisted PCAC 
quark mass mpcAC- The small discrepancy observed at j3 = 0.67 (left panel) between extrapolations from 
positive and negative quark masses is probably a small effect of the first order phase transition. For j3 — 0.74 
(right panel), extrapolations from both sides give consistent results. An alternative way to fix the critical 
mass is also plotted (see for details). 



Indeed using an off diagonal splitting, where the degenerate quark doublet has a different flavour 
orientation from the splitting between the quarks, in [pO[] it was shown that the determinant for non 



degenerate quarks is real and positive (see [47] for alternative formulations of tmQCD including 
a non degenerate doublet). The only restriction of the construction in [^] is on the value of 
the ratio between renormalization constants of the pseudoscalar and scalar current. To give an 
example, fixing the values of the renormalized strange and charm quark masses, gives the following 
constraints 

~ l.SGeV ~ O.lGeV 0.875. (3.4) 

At this conference first results with dynamical Nf = 2 + 1 + 1 twisted quarks have been presented 
[76]. The simulations of the N{ = 2 + 1 + 1 theory are performed by a polynomial hybrid Monte 
Carlo algorithm (PHMC)[81 ]. The structure of the algorithm goes along the lines indicated in [^]. 
At this conference another variant of the PHMC to include the two non degenerate twisted quarks 
has been presented [83]. 



4. Further results 

In this section I summarize further results concerning tmQCD. In [ ^4| ] it has been presented 
a strategy to compute Bk and matrix elements related to the A/ = 1/2 rule without mixing with 
operators with wrong chiralities, retaining all the properties of automatic O(fl') improvement. The 
strategy is based on the usage of a mixed action (OS for valence quarks and tmQCD for sea quarks) 
[pT]]. In the quenched approximation Bk has been computed [ p5[ ] in the continuum limit, us- 
ing clover improved tmQCD and a non-perturbative renormalization without mixing (in the SF 
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scheme). A strategy to compute with tmQCD, along the Unes of pq], has been proposed in 
^T\\, and along the lines of Q] in [^]. In [47] a strategy, based on clover improved tmQCD with 



A'f = 4, has been proposed to compute the renormalization of ^ ^ tt matrix elements. In [89] the 
effect of a twisted mass term of the low-lying modes of the Wilson-Dirac operator and a remnant 
of the index theorem for twisted mass fermions has been discussed. 



5. Conclusions 

Several lessons come from quenched studies of tmQCD. With a particular and field theoreti- 
cally well founded definition of the critical mass, automatic 0(a) improvement is effective till small 
pion masses (nijc = 272 MeV), and the residual O(a^) cutoff effects are small. The bending phe- 
nomenon just results from big cutoff effects, that are reproducible with ;^PT at finite lattice spacing. 
The bending phenomenon is not present even at finite lattice spacing with a suitable choice of the 
critical mass. The flavour breaking is an issue and it has to be investigated with dynamical simu- 
lations. We have indications of the existence of an Aoki phase for quenched Wilson fermions at 
lattice spacings around a ~ 0.1 fm. 

To perform dynamical simulations at small pion masses, algorithmic improvements are crucial, 
and now new algorithms allow to have efficient and performant simulations with Wilson fermions 
and most probably with staggered fermions. 

We have a much better understanding of the phase structure of dynamical Wilson fermions. 
A theoretically well founded action (tlSym gauge and tmQCD fermion action) allows to perform 
dynamical simulations with Nf = 2 at pion masses smaller then 300 MeV starting from a lattice 
spacing a ~ 0.1 fm, allowing matching with ;^PT, and simulation with A'^f = 2 + 1 + 1 flavours are 
just starting. 

In many cases it has been shown that the renormalization properties of local operators related 
to very important phenomenological quantities, is continuum like. 

Although presently not all aspects of tmQCD are fully investigated, tmQCD is an attractive and 
powerful discretization of lattice QCD, and it certainly belongs to the pool of well founded fermion 
actions that ought to be used to control the continuum hmit of physical quantities of interest. 
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